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ABSTRACT 


In  order  to  maintain  the  amplitude  of  mechanical  vibrations 
in  a  solid  constant,  it  is  necessary  to  supply  energy  to  the  system. 
This  dissipation  of  energy  may  occur  via  many  distinct  mechanisms. 

The  present  work  was  an  experimental  study  carried  out  in 
order  to  determine  the  mechanism  which  causes  the  anomalously  large 
dissipation  observed  in  iron,  at  low  temperatures  for  vibration 
frequencies  in  the  KH^  range,  by  Bruner  (i960),  Heller  (l9'6l)  and 

4O1 

Verdini  (1962).  It  will  be  shown  that  it  is  possible  to  explain 

this  behavior  as  being  due  to  the  damping  of  dislocations  by  the 

conduction  electrons  in  the  metal,  assuming  a  dislocation  density 
7  ..o 

of  10'  cm  ^  and  an  average  length  for  these  dislocations  of 
4.3  x  10"*  cm.  These  values  ar*  in  good  agreement  with  those 
found  by  other  techniques. 

In  the  course  of  this  work,  the  damping  of  vibrations 
was  measured  in  two  other  ferromagnetic  metals,  cobalt  and  nickel. 

The  results  for  cobalt  showed  evidence  of  the  unstable  hexagonal 
phase  of  this  metal  which  has  previously  been  seen  by  other  methods 
(Bozorth,  1951) 
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CHAPTER  I 
INTRODUCTION 


In  an  ideal  elastic  solid  Hooke’s  law  is  obeyed,  that  is 
each  stress  component  is  a  linear  real  function  of  the  strain  compo¬ 
nents  alone.  Within  this  law,  there  is  no  allowance  for  damping  of 
vibrations  or  traveling  waves.  In  a  real  solid,  however,  vibrations 
and  x^aves  are  attenuated  and  the  relation  between  stress  and  strain 
components,  although  linear,  must  contain  terms  in  stress  and  strain 
rate  (Zener,  1948),  This  non  ideal  behavior  stems  from  the  fact  that 
phonons,  the  quanta  of  vibrational  energy,  interact  with  electrons, 
dislocations,  point  defects  etc.,  in  the  solid,  which  absorb  energy  in 
the  process  giving  rise  to  attenuation.  Vibrations  in  a  solid  can 
also  set  up  thermal  currents  and  since  this  heat  flow  is  resistive, 
elastic  energy  is  absorbed  and  this  produces  still  another  form  of 
attenuation.  Heat  plays  another  role  in  attenuation  in  that  vibra¬ 


tional  phonons  can  be  scattered 
be  scattered  by  grain  boundaries 
the  scattering  processes  produce 
or  less  limited  to  the  megacycle 


by  thermal  phonons.  Phonons  can  also 
in  polycrystaline  solids.  Each  of 
s  attenuation,  however  they  are  more 
frequency  range. 


Each  of  the  attenuation  processes  mentioned  above  affords  a 


method  of  experimental  study  on  the  parameters  involved.  For  instance 
it  is  possible  to  measure  the  diffusion  constant  of  point  defects  such 
as  vacancies,  interstitial  or  substitutional  impurities.  A  study  of 
damping  due  to  dislocations  can  lead  to  a  value  for  the  height  of  the 


Peierls  barrier  or  the  binding  energy  of  impurities  to  dislocations. 
Measurements  of  the  electron  phonon  interaction  will  yield  the  shape 
of  the  Fermi  surface  of  a  metal  or  the  energy  gap  of  a  superconductor. 
Comprehensive  reviews  of  the  various  theories  involved  and  their  rela¬ 
tion  to  experiment  have  been  given  by  Ilason  (1966)  and  Bhatia  (1967). 
The  method  of  study  usually  involves  the  measurement  of  the  temperatur 
and  frequency  dependence  of  the  vibration  damping  and  of  the  elastic 
modulus . 


The  present  work  was  a  study  of  the  anelastic  behavior  of 
poly cryst aline  iron  with  emphasis  021  the  lot-;  temperature  region.  The 
method  used  was  to  measure  the  damping  of  free  vibrations  set  up  in 
circular  plates  as  a  function  of  several  parameters,  the  most  direct 
of  which  was  temperature.  The  purpose  of  the  experiment  was  to  inves¬ 


tigate  and  explain  the  large  increase  in  attenuation  at  low  tempera 
tures  observed  by  Brunei*  (i960),  Heller  (l$)6l)  and  Verdini  (1962). 
Measurements  were  carried  out  over  the  temperature  range  from  4,2 0 


ft 


3 


to  room  temperature  at  frequencies  ranging  from  10  KH  to  170  KH^. 
The  dependence  of  the  anomalous  damping  on  magnetic  field,  impurity 
content,  and  thermal  and  mechanical  treatment  was  also  studied.  The 
dependence  of  damping  on  temperature  in  cobalt  and  nickel  was  also 
measured  so  that  the  behavior  in  iron  could  be  compared  with  that  of 
other  ferromagnetic  metals. 
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CHAPTER  II 


BASIC  THEORY  OF  AK5LAST IC IIY 


A*  Ideal  Elasticity 

The  formulation  of  the  theory  of  elasticity  began  with 
Robert  Hooke’s  law  which  reads: 

F  =  KD  ,  (2.1) 

- ^  ^ 

where  F  is  a  force  acting  on  a  body  to  produce  a  deformation  D,  The 

proportionality  constant  M  is  called  the  elastic  modulus  of  the  body. 

The  formulation  however  required  that  the  force  be  applied  slowly 

enough  that  the  whole  of  the  body  be  aware  of  its  presence.  This 

limitation  stems  from  the  fact  that  a  mechanical  pulse  has  a  finite 

velocity  namely,  the  speed  of  sound  in  the  material  under  consideration. 

Another  fact  not  taken  into  account  by  Hooke’s  law  is  that  a  force 

suddenly  applied  to  a  body  may  set  up  vibrations  in  this  body.  In  order 

for  this  law  to  be  valid  then,  the  force  must  be  applied  over  a  time 

which  is  longer  than  the  period,  of  the  lowest  frequency  normal  mode  of 

the  specimen. 

Voigt  (1910)  overcame  this  difficulty  by  generalizing  Hooke's 


- - ------  - - - — i - n« 
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law  to  consider  all  forces  including  the  rapidly  varying  ones.  To  do 
this  one  divides  the  material  up  into  small  rectangular  parallelepiped 
volumes  as  illustrated  in  figure  2.1,  of  dimensions  dx^  dx-> ,  and  dx^ 
and  considers  only  the  forces  and  deformations  in  these  volumes  due  to 
the  surrounding  material.  In  the  limiting  case  of  dx_.  tending  to  zero, 
the  period  of  the  lowest  normal  mode  will  tend  to  zero  and  any  force 
applied  to  the  elemental  volume  will  be  instantly  felt  in  the  form  of 
a  deformation.  The  components  of  the  force  per  unit  area  on  each  face 
of  this  volume  are  called  the  stress  components  and  are  noted  cx^.  The 
indices  indicate  the  face  on  which  the  force  acts  and  its  component 
respectively.  We  now  define  the  strain  components  ^  as  •; 


£ ij  =  |[(3si/axj)  -i-  (asyaxi)]  , 


whore  s^(x)  are  the  components  of  the  displacement  of  particles  at  x. 

law 

Both  a.  .  and  £ .  .  behave  as  second  rank  tensors  and  Hooke’s  can  be  writ- 
-*-<3  •*- J 

ten  as  ; 


ai j  ~  cijkl“kl 


(2.2) 


where  the  coefficients  e-  1  are  called  the  elastic  constants  of  the 

J- JK-L 


medium. 


I 
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Figure  2.1 


Rectangular  parallelepiped  portion  of  a  medium  with  stress  components 
illustrated,  on  one  face  (after  Bhatia  1967). 


■gyggc*-  vrga-yag-.  sj«3re  tts-.  ^a^rwutrffcar  •>  v.t  ri^'nts^vme^r^gaw 
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With  equation  (2.2)  in  mind,  the  equation  of  motion  of  the 
particles  at  a  point  can  be  written  as : 


as 


.2 

a  s 


at2  ijk-L  c>Xj5x^ 


(2.3) 


where  s^  is  again  the  component  of  displacement.  This  is  a  wave  equa¬ 
tion  and  has  the  usual  plane  and  standing  wave  solutions  and  these 
waves  have  undiminishing  amplitudes  in  time  or  distance  travelled. 


B,  Real  Solid s 


(1)  Non  Elastic  Behavior 

The  solid,  described,  by  equation  (2.2)  behaves  mechanically 
as  shown  in  figure  2.2(a)  that  is  a  force  instantaneously  applied  at 
a  time  t^  produces  an  instantaneous  deformation  which  remains  constant 
as  long  as  the  force  does,  .The  deformation  will  then  go  instantaneously 
to  zero  when  the  force  does. 


all ,  the 
limit , 


A  real  solid  behaves  in  a  very  different  manner, 
magnitude  of  the  force  is  limited  to  what  is  called 
Any. force  greater  than  this  will  produce  a  permanent 


First  of 
the  elastic 
deformation 


in  the  solid,  a  phenomenon  called  plasticity  of  which  no  account  is 
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Figure  2.2 


The  mechanical  behavior  of  (a)  the  ideal  elastic  solid  and  (b)  the 
standard  linear  solid  (Zener,  1948). 


t.z  time 


taken  in  Hooke’s  law.  Even  -while  still  within  the  elastic  limit,  a 
real  solid  does  not  behave  elastically  but  rather  in  the  manner  de¬ 
picted  in  figure  2.2(b). 
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(2 )  An elastic it v 


In  order  to  describe  the  behavior  within  the  elastic  limit, 
one  must  modify  equation  (2.2)  so  that  it  includes  terms  in  the  var¬ 
ious  orders  of  the  time  derivitive  of  the  stress  and  strain  components 
(Poynting  and  Thomson,  1905 )•  A  very  good  approximation  to  this  is 
given  by  the  equation  containing  only  the  first  time  derivitives  of 
stress  and  strain.  The  one  dimensional  equation  for  this  case  is  writ¬ 
ten: 


a-^a  •+  a^cj  — 


+  b26 


(2.4) 


This  equation  contains  three  independent  constants  which  are  chosen 
such  that  the  equation  can  be  written, 


o  +  ry  =  1 ig(t  +rje). 


(2.5) 


Solids  which  obey  this  equation  are  called  standard  linear  solids. 
We  now  suppose  that  a  stress  cr 


is  suddenly  applied  at  t  = 0 


■ 


. 
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and  is  held  constant  thereafter;  then  equation  (2.5)  becomes; 


cr 


+ 


t  >  0, 


which  when  rewritten  is, 


€  +  e-:1  =  (%^)-1a0  . 


The  solution  to  this  equation  is, 


£(t)  =  Mfc,  +  (SD  -  , 


(2.6) 


showing  that  the  behavior  of  the  standard  linear  solid  is  as  that 
shown  in  figure  2.2(b)  for  a  real  solid.  It  is  seen  that  the  strain 
c  is  equal  to  at  a  time  t  equal  to  zero  and  then  relaxes  to  its 
maximum  value  il  as  t  becomes  large.  Since,  after  all  relaxation 
has  oc cured, 


Lira  a 

t-»  OO  £ 


(2.7) 


the  constant  lh,  trill  be  called  the  relaxed  modulus. 
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tv  W9 


In  order  to  find  the  physical  significance  of  X  and 
suppose  a  case  where  €  and  6  are  identically  zero.  Then  equation 
(2.5)  reads, 


a  +  Xa  =  0, 


whose  -solution  is 


a(t)  —  a(0)e~^x 


This  indicates  that  is  the  relaxation  time  of  the  stress  under 
conditions  of  constant  strain;  similarly  is  the  relaxation  time  of 
the  strain  under  constant  stress. 

Now  if  we  consider  a  variation  in  stress  ha  inducing  a  varia. 
tion  in  strain  Ac  over  a  vanishingly  small  time  6t,  we  find  from  equa¬ 
tion  (2.5)  * 

A  a  =  Mr  «  (2.8) 


The  quantity  M  , t”0/ will  be  note  11^  and  is  called  the '  unrelaxed  mqdu- 
lus  since  it  relates  changes  in  stress  and  strain  which  occur  in  such  a 
short  period  of  tine  that  no  relaxation  has  time  to  occur. 


Since  the  present  investigation  is  concerned  with  the  behavior 
of  solids  under  periodic  stress  we  must  investigate  the  properties  of 
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equation  (2.5)  with  this  in  mind.  In  substituting  into  this  equation 
the  following  solutions: 


a(t)  ~  aoe 


iox 

> 


and 


fe'(t) 


f 

“O 


♦ 


it  is  found  that 


o 


o 


=  1)1G 

O 
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where  >W  is  complex  and  given  by, 


In 


1  +  iw  ^ 

_ _  a 

1 


(2.9) 


The  fact  that  the  modulus  relating  a  and  6  is  complex  indi¬ 
cates  that  the  strain  lags  behind  the  stress  by  a  phase  which  we  shall 
denote  6.  It  is  easily  seen  that 


tan 


Ira  Tfl 
Re  >77  ? 


which  is  evaluated  from 


(2.9)  to  obtain: 


tan  6  = 


6>(ta-r£) 

i+«2(Vf) 


3 


(2.10) 


13 


If  mechanical  vibrations  are  set  up  in  a  specimen  by  a  periodic  stress; 


a  =  a  cos  <ut, 
o 


the  resulting  strain  will  be,  in  general, 


6  =  60  cos  (a)t  -  6). 


The  rate  at  which  energy  will  be  dissipated  per  unit  volume  during  such 
vibrations  is  the  mean  value  of  j  ,  then  the  total  power  dissipation  for 
a  sample  will  be: 


P  = 


dv, 


sin 


o 


dv 

o  o 


(2.11) 


If  tan  6  is  small  then  the  vibrational  energy  per  unit  volume  can  be 

taken  as  4a  £  so  that  the  total  energy  is 
A  o  o 


b 


JL 

2 


o  c  dv, 
o  o 


v/ 


(2.12) 


A  standard 


quantity 


rr 


will  now  be  defined  in  the  usual  way  for  a 
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driven  oscillator: 

q-1  _  Power  dissipation 

w  X  energy  of  vibration  * 

where  cj  is  the  angular  frequency  of  the  vibrations.  For  the  case  of 
the  standard  linear  solid,  it  is  seen  from  equation  (2 .11)  and  (2.12) 
that 

Q"1  =  sin  (2.13) 

and  since  tan  5  was  assumed  to  be  small, 

Q-1  «  tan  6  .  (2.14) 

Tan  6  has  been  evaluated,  in  terms  of  the  relaxation  times  in  equation 

(2.10). 

Equation  (2.14)  states  that  the  normalized  power  dissipation 
in  a  standard  linear  solid  is  non  zero  since  tan  6  is  non  zero,  that 
is  vibrations  are  damped  due  to  the  fact  that  some  of  their  energy  is 
being  dissipated  to  the  medium.  Since  a  real  solid’s  behavior  is 
approximated  by  the  standard  linear  solid,  vibrations  should  also  be 


' 


. 
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damped  in  their  case,  and  it  is  a  well  known  experimental  fact  that 
this  is  true. 

It  has  been  seen  that  the  value  of  is  a  measure  of  the 
amount  of  damping  that  vibrations  in  a  real  solid,  are  submitted  to. 
Experimentally,  Cf1  could  be  measured  by  measuring  the  phase  lag  of 
the  induced  periodic  strain  relative  to  the  stress  applied.  This  how¬ 
ever  is  a  difficult  task  since  the  angle  5  is  usually  a  very  small 
quantity.  In  order  to  avoid  this,  Q“^  must  be  related  to  other  quan¬ 
tities  more  easily  measured  in  the  laboratory. 

However  before  discussing  the  measurement  of  Q”^*  we  shall 
return  to  equation  (2.10),  which  can  now  be  written: 


&  tan  6  = 


w(t  -  TJ 

a  t 

i  +</(Var£  j 


(2.3.5) 


to  see  some  of  its  other  properties,  lie  first  define  the  geometric 
mean  of  the  two  relaxation  times  and  T'  as, 


?  -  c 


T  ) 

o' 


1/2 


j 


(2.16) 
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and  that  of  the  two  elastic  moduli,  M  and  M0  as, 

*  u  R  * 


M  =  ( 


1/2 


(2.17) 


Using  these  definitions,  equation  (2,15)  Diay  be  rewritten  as, 


Q 


M  -Ma 
u  1 


CO  L 


(CO) 


'ti  1  +  (ij  r )< 


(2.18) 


which  gives  the  frequency  dependence  of  Q~  ,  Equation  (2,18)  is  de- 
picted  in  figure  2.3.  The  maximum  value  of  Q  which  occurs  at  U)v  =1, 
is  given  by: 

•1 


Q. 


M  -1L 
u  R 


max 


2M 


It  was  shown  in  equation  (2.9),  for  the  case  of  periodic 
stress  the  elastic  modulus  is  complex.  It  is  of  value  to  define  a  real 
modulus  M  which  can  be  measured  experimentally.  We  choose  to  define 


s.t  as, 


M 


co 


[ro  /tT  j' 


-1 


that  is  TI,,  is  the  ratio  of  the 

£ j ) 


stress  to  that  part  of  the  strain  which 


Figure  2.3 


The  frequency  dependence  of  Q' 


and  M^.  (after  Zener,  1948) 
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is  in  phase  with  it  (Zener,  1948,  pp . 4? ) .  From  equation  (2.9)  it  is 
seen  that 


.  2^2 
1  +  0  1  a 

-i  2  ~  ~ 

1  +  lT  l  l 
O  t 


9 


Then  by  using  equation  (2.8)  the  final  result  is: 


\  ^ 

l  + 


I 


which  has  the  limiting  values, 

;  cJt  «  l 

—  a 

;  ut  »  1 

This  behavior  is  shown  in  figure  2.3. 

The  theory  presented  above  was  developed  for  the  standard 
linear  solid.  It  covers  damping  of  vibrations  characterized  by  a 
single  relaxation  time  X .  It  is  a  very  good  approximation  to  the 
behavior  of  red  solids  in  which  the  damping  is  governed  by  one  relaxa- 


, 
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tion  time.  It  is  possible  as  well  to  develop  a  theory  which  covers 
damping  with  a  spectrum  of  relaxation  times  by  superposition  of  the 
results  of  each  relaxation  time  in  the  spectrum.  This,  however,  does 
not  give  a  great  deal  of  insight  into  the  physical  phenomena  involved 
and  for  this  reason  most  data  is  today  analysed  from  first  principles 
rather  than  the  relaxation  spectrum  of  a  phenomenological  theory.  For 
example,  in  the  theoretical  interpretation  of  damping  due  to  electron 
phonon  interaction,  one  considers  the  rate  of  energy  dissipation  due 
to  the  disturbance  of  the  viscous  electron  gas  by  the  vibration  phonons 
(Bhatia,  1967). 


(3)  Definition  of  Parameters  Related  to  Elastic  Energy 


Dissipation 


^-1 


It  is  now  necessary  to  return  to  the  parameter  Q  .  From 
the  defining  equation  it  is  known  that: 


Q-X  = 


Power  dissipation 
u)  x  energy  of  vibration 


(2.19) 


If  we  define  ZiE  as  the  energy  dissipated  per  cycle  of  vibration  of 
total  energy  E,  then, 
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A2  __  2tt  Power  dissipation 

E  ~ 


(2.20) 


On  comparing  equations  (2.19)  and  (2.20)  it  is  found  that, 


1  /I  AS 

~  TT  2  E 


); 


(2.21) 


thus  Q  ^  has  been  related  to  the  energy  dissipated  during  one  half 
cycle  of  vibration.  The  energy  S  must  now  be  related  to  the  amplitude 
of  vibration  which  is  a  measurable  quantity.  Since  elastic  energy  is 
dissipated,  the  amplitude  of  free  vibrations  must  decrease  in  time. 
Suppose  the  amplitude  at  the  beginning  and  the  end  of  a  period  of  time, 
equal  to  the  period  of  the  vibration,  is  known.  These  amplitudes  will 
be  noted  A  and  A-AA  respectively.  A  new  parameter  called  the  logarith¬ 
mic  decrement  is  now  defined  as, 


logarithmic  decrement  =  In 


(2.22) 


Since  the  elastic  energy  of  the  vibration  is  proportional  to  the  square 


then, 


A 


,1/2 


A-AA  (E_lE)l/2 


(2.23) 
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where  AS  was  defined  for  equation  (2,21),  If  AE  is  small  it  is  easily- 
shown  that, 


log,  dec.  -  In 


AE 


iUM  2  E 


(2.24) 


In  comparing  equations  (2.21)  and  (2.24)  it  is  seen  that  ! 


log.  dec  =  ttQ”\ 


(2.2 5) 


From  equation  (2. If)  it  is  seen  that  is  dependent  on  frequency  and 
the  two  relaxation  times  f  and  T^.  Since  the  frequency  is  independent 
of  amplitude  and  it  is  assumed  that  the  relaxation  times  are  also  inde¬ 
pendent  of  amplitude  then  equation  (2.2f)  can  be  written  in  the  form$ 

l-l 

A(t)  =  A(0)  exp ,  ( -  )wt  o  (2.2 6) 


This  shows  that  in  a  standard  linear  solid  the  amplitude  of  free  vibra¬ 
tions  decays  exponentially.  It  is  a  well  established  experimental  fact 
that  in  real  solids  for  the  most  part  the  decay  of  vibrations  is  expo¬ 
nential.  This  suggests  a  method  for  measuring  Q”~;  that  is  by  measur¬ 
ing  A(t)  for  free  decay  and  evaluating  the  exponential.  This  method 
is  practical  for  value  of  the  dissipation  which  give  decay  times  that 
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are  long  enough  to  be  measured  accurately. 

Suppose  we  have  a  sample  which  is  in  forced  vibration  near 
its  resonant  frequency.  It  is  possible,  in  this  case,  to  evaluate 
the  dissipation  from  the  shape  and  frequency  of  the  resonance  peak. 

It  is  assumed  that  the  sample  can  be  represented  by  a  system  with  one 
degree  of  freedom  and  with  an  anelastic  restoring  force,  f,  similar  to 
the  standard  linear  solid.  If  the  displacement  from  equilibrium  of 
this  system  is  noted  x,  then  the  equation  of  motion  can  be  written, 


px  +  f 


V 


±0)t 
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where  p  is  the  inertia  of  our  system  and  Fo  the  amplitude  of  the  driv¬ 
ing  force.  The  restoring  force  is  related  to  the  displacement  by: 

f  =  >n  x  =  M(l  -ri  tan  o)x  > 


where  'M  is  the  complex  elastic  modulus.  Assuming  a  solution  of  the 


form: 


x 


itot 


xoe 


9 


and  substituting  into  the  equation  of  motion,  it  is  seen  that, 
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/M 


x.  = 


(n  P'JS 


(2.2?) 


+  i  tan 


The  amplitude  of  the  displacement  is  seen  so  have 


ci.  j.ucOs.J 


ui  at  the 


frequency  <r  given  by: 


“o  =  (f ' ' 


ana  to  have  half  this  maicimrwi  value  when 


cJ2 

o 


+  „  /  d 

_  V  J 


m  o 


From  equation  (2.1 h)  it  is  knovm  that 


—  o  an  o  . 


-1 

therefore,  Q  is  related  to  the  width  of  the  resonance  peak  at  half 
maximum,  denoted  by  AdO,  in  the  following  way: 


Q' 


A  6) 

' 


(2.28) 


CHAPTER  III 


EXPERT!  EHTAL  TECHNIQUE 

A*  Samples  and  mounting 

The  material  to  be  investigated  in  this  work  was  polycrysta- 
line  iron,  more  specifically  Joins on,  Matthey  spectrographically 
standardized  iron.  It  was  purchased  in  the  form  of  a  plate  of  3  mm 
thickness  from  which  disks  were  cut  with  diameters  of  ^0 '  and  mm. 

The  resonance  technique  was  used  to  investigate  the  acoustic  damping  in 
this  material;  that  is  the  samples  were  excited  to  vibration  in  their 
flexural  normal  modes  and  depending  on  circumstances  (either  the  rate 
of  decay  of  free  vibration  or  the  width  of  the  resonance  peak  was  meas¬ 
ured.  The  method  of  mounting  the  samples  in  the  cryostat  was  that  de¬ 
scribed  by  fordoni,  Huovo,  and  Verdini  (1959). 

The  point  behind  this  method  is  minimizing  any  damping  of 
vibrations  in  the  specimen  due  to  the  surroundings.  One  cause  of  exter¬ 
nal  damping  is  acoustic  radiation  from  the  sample  to  the  surrounding 
air.  This  is  effectively  nullified  by  placing  the  sample  in  an  evacu¬ 
ated  chamber.  Jill  the  measurements  presented  here  were  done  at  pres- 
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sures  of  less  than  10”^  mu  Hg,  The  next  problem  arises  with  supporting 
the  sample  in  the  chamber.  A  schematic  diagram  of  the  mount  used  is 
shown  in  figure  2.1.  It  consists  of  resting  the  sample  on  three  sharp 
pins  which  come  into  contact  with  it  on  a  nodal  circle  and  are  separated 
by  120°  on  that  circle.  In  order  to  prevent  the  sample  from  moving  on 
the  pins,  three  shallow  holes  were  drilled,  at  the  points  of  contact  of 
the  pins.  The  fact  that  the  sample  was  contacted  only  at  three  sharp 
points  located  on  a  nodal  circle  minimized  the  amount  of  elastic  energy 
transferred,  from  the  sample  to  its  mount. 

These  pins  served  other  functions ,  as  well  as  being  a  low 
damping  support.  One  pin  was  steel  and  provided  a  ground  connection 
for  the  sample  to  the  cryostat  as  illustrated  by  "3"  in  figures  3.1  and 
3.2.  The  other  two  pins  (1  and  2  in  figures  3*1  and  3«2)  were  in  the 
form  of  fine  wires  of  about  0.33  mm  in  diameter,  insulated  from  ground 
and  were  used  as  a  thermocouple  with  the  sample  forming  the  junction. 

The  material  used  for  the  thermocouple  was  gold  -b  2,1  At  Co  refer¬ 
enced  to  copper.  This  was  chosen  because  it  is  the  most  sensitive  low 
temperature  thermocouple  readily  available  ( Powell  et  al. ,  1961). 

B,  The  Excitation  and  Detection  of  Vibrations 

The  vibrations  in  the  sample  were  excited  by  an  electrostatic 


Figure  3.1 


Schematic  diagram  of  the  sample  mount. 


C  -  nodal  circle 
E  -  driving  electrode 
P  -  support  pins 

1  and  2  -  thermocouple 

3  - 

S  -  sample 


ground  pin 


sn 


Figure  3*2 


Schematic  diagram  of  the  Cryostat  showing  the  sample  mounted 
(see  figure  for  labels) 

A,  B  and  C  -  level  adjustments  for  support  pins 
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force  from  an  electrode,  in  the  form  of  a  plate  of  about  6  mm  in  diam¬ 
eter,  placed  directly  above  and  parallel  to  the  surface  of  the  sample 
at  a  distance  of  a  few  tenths  of  a  millimeter.  This  is  illustrated  by 
E  in  figure  3«1  and  j.2a  The  distance  of  the  electrode  from  the  sample 
was  controlled  by  turning  it  in  a  finely  threaded  mount  in  the-  sample 
chamber.  The  turning  was  accomplished  via  a  small  diameter,  thin  wall, 
stainless  steal  tube  connected  directly  to  the  electrode  and  exiting 
the  cryostat  up  the  main  pumping  line  through  a  rubber  vacuum  seal. 

The  control  rod  was  insulated  from  ground  and  served  as  an 
electrical  feed-through  to  the  electrode.  The  excitation  and  detection 
circuits  connected  to  this  rod  are  shown  in  the  block  diagram  of  fig¬ 
ure  3*3 •  In  order  to  excite  the  vibrations  in  the  sample,  it  was  kept 
at  ground  potential  by  means  of  pin  3  in  figure  3.2  and  a  sinosoidal 
voltage  was  applied  to  the  electrode  through  the  control  rod.  In  order 
to  calculate  the  force  on  the  sample  due  to  a  voltage  V  applied  to  the 
electrode  it  is  necessary  to  consider  the  system  as  a  parallel  plate 
capacitor.  If  the  capacitance  of  the  system  is  C  and  the  charge  on  the 
plates  is  given  by  q,  then  the  energy  W  stored  in  the  system  will  be 
given  by: 


W 


1 

C 


0 


q’dq>  , 


Figure  3.3 


Block  diagram  of  the  experimental  set 
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therefore , 


r2 

W  =  i  Si 
2  C 


1 


Since,  for  a  parallel  plate  capacitor,  C  is  given  by: 


C  = 


9 


where  (£q  is  the  permittivity,  A  the  area  of  the  plates  and  d  the  dis¬ 
tance  that  separates  them,  then  the  force  1  on  the  sample,  using* 


F 

will  be  given  by- 


and  C  =  ^  , 
V 


F 


(3.1) 


If  the  applied  voltage  V  is  of  the  form* 


V  = 


T  iot 


9 


then,  since  the  resulting  force  F  is  proportional  to  the  square  of  this 
voltage,  it  will  be  given  by: 
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F  =  F  e' 


2  lot 


which  is  an  oscillating  force  of  frequency  2^.  This  force  will  then 
set  up  flexural  vibration  in  the  sample  of  this  same  frequency  20. 

The  driving  voltage  came  from  an  oscillator-amplifier  com¬ 
bination  whose  output  level  could  be  varied  from  0  to  250  volts  over  a 

3oo 

frequency  range  of  2  KH  to  KH  .  The  voltage  amplifier  x-xas  built  in 
this  laboratory  and  x-ias  similar  to  that  described,  by  IJuovo  (Ifol) . 

The  oscillator  was  constructed  to  have  a  x«ride  band  and  high  stability 
according  to  the  circuit  shown  in  figure  3*^«  The  tuning  x-ras  such  to 
allow  a  change  of  the  output  frequency  by  less  than  0.1  H  over  the 
entire  frequency  range .  The  frequency  stability  was  f ound  to  be  better 
than  0.1H  after  a  few  hours  of  operation.  The  fine  tuning  ability  was 
achieved  by  using  a  100  to  1  mechanical  reduction  on  the  controls  of 
the  variable  capacitors  A  and  E.  In  summary ,  the  vibrational  state  of 
the  sample  was  controlled  by  the  frequency  and  amplitude  of  the  driving 
voltage  and  the  distance  of  the  electrode  from  the  specimen. 

The  fact  that  the  electrode  and  the  sample  formed  a  parallel 
plate  capacitor  was  taken  advantage  of  in  order  to  measure  the  amplitude 
of  the  excited  vibrations.  Y/hen  the  sample  was  vibrating  its  average 
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Figure  3.^ 


Cathode-coupled  oscillator  used  for  the  driving  voltage 
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distance  from  the  electrode  varied  inUsinosoidal  way  thereby  causing 
a  variation,  C,  in  the  capacitance  of  the  system.  Considering  that: 


where  d  is  the  electrode  to  sample  distance  then  a  sinoso.idal  variation 
in  d  of  amplitude  Ad  will  produce  a  corresponding  variation  in  C  whose 
amplitude  will  be  proportional  to  Ad.  The  experimental  problem  is  now 
to  measure  AC  in  order  to  get  a  relative  measure  of  the  vibrational 
amplitude.  In  referring  to  chapter  II,  section  B(3)  it  is  seen  that  a 
measure  of  relative  amplitude  suffices  to  specify  the  magnitude  of 
The  measurement  of  AC  is  done  by  including  the  sample-electrode  capac¬ 
itor,  C,  in  an  antiresonant  discriminator  circuit  whose  input  is  con¬ 
nected  to  an  rf.  oscillator  with  a  10  Mri^  output.  Because  of  the  varia¬ 
tion  in  C  the  ''resonant'1  frequency  of  the  discriminator  is  modulated 
causing  an  amplitude  modulation  of  the  10  MH  signal  across  the  dis- 
criminator.  The  frequency  of  the  side  bands  is  the  same  as  that  of  AC 
and  their  amplitude  is  proportional  to  it,  thus  to  the  amplitude  of 
vibration.  This  circuit  was  described  in  a  paper  by  Kuovo  (196l)  along 
with  the  demodulator  which  produced  an  audio  frequency  signal  corre¬ 
sponding  to  the  sidebands  of  the  modulated  rf.  signal,  and  therefore 
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the  vibration  of  the  specimen. 


This  signal  was  then  put  through  a  variable  high-pass  filter, 
set  at  about  5  KH  below  the  frequency  of  vibration,  to  cut  out  low 
frequency  noise  such  as  room  vibrations  and  any  leak-through  which 
might  have  occured  from  the  driving  signal.  The  signal  was  then  routed 
through  a  tuned  amplifier  system  whose  circuit  is  illustrated  in  fig¬ 
ure  3.5.  The  tuned  stage  of  this  amplifier  has  a  Q  which  is  much  lower 
than  the  mechanical  one  in  order  to  avoid  interference  with  the  measure 
ments.  -The  amplified  signal  was  then  fed  into  a  logarithmic  level  re¬ 
corder  (Bruel  and  Kjaer ,  model  2305)  and  vacuum  tubs  volt  meter.  This 
signal  was  also  used  as  the  y  input  of  an  oscilloscope  whose  x  input 
was  connected  to  the  driving  signal.  The  2  to  1  Lissajou  figure  dis¬ 
played  served  as  a  monitor  for  the  detected  and  driving  signals.  The 
driving  signal  frequency  was  measured  to  an  accuracy  of  +0.1  Ha  by  a 
frequency  counter  connected  to  the  output  of  the  voltage  amplifier. 

C.  Temperature  Measurement 


After  the  frequency  and  relative  amplitude  of  the  vibrations 
are  known,  a  third  parameter  is  needed,  namely  the  temperature  of  the 
specimen.  The  thermocouple  described  earlier  provides  this  information 
The ’reference  junction  was  placed  in  an  ice  bath  whose  temperature  was 
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Figure  3.5 


Preamplif ier  and  tuned-amplifier  used  in  the  detection  eireui 
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monitored  by  a  mercury  thermometer  with  a  resolution  of  0.1°  C.  The 
calibration  of  this  thermocouple  was  derived  from  data  given  by  Powell, 
Bunch,  and  Corruchini  (I96I)  for  a  gold  plus  2.1  At.  $  cobalt  versus 
copper  thermocouple.  To  do  this  the  thermocouple  e.m.f.  was  measured 
at  the  boiling  point  of  oxygen  and  helium  and  the  ratio  between  these 
values  and.  those  given  by  the  above  authors  were  compared  and  found  to 
be  equal  to  within  0,1$  so  it  was  assumed  that  it  would  be  a  constant 
throughout  the  temperature  range.  From  this  assumption  a  table  was 
prepared  giving  thermocouple  e.m.f.  versus  temperature  at  intervals  of 
0.5°  K  from  1  to  3^0°  K.  The  output  of  the  thermocouple  was  measured 
by  standard  potentiometric  means  and  compared  to  this  table  by  inter¬ 
polation  to  arrive  at  the  temperature  of  the  sample.  Since  the  whole 
of  the  sample  formed  the  junction  of  the  thermocouple,  the  temperature 
measured  in  this  way  was  the  average  temperature  of  the  sample.  The 
precision  of  this  measurement  was  estimated  to  be  +  0.2°  K  from  room 
temperature  to  about  15°  K.  Below  this  temperature,  the  error  in¬ 
creased  to  its  maximum  value  of  ±  0.5°  K  at  helium  temperatures. 

D.  Operating  Procedure 

In  order  to  demonstrate  the  use  of  the  equipment  described 
above  it  is  instructive  to  go  through  the  operations  of  a  typical  data 
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run.  The  first  step  is  to  set  the  sample  on  the  support  pins  and  to 
adjust  the  pin  holders  (A,B  and.  C  in  figure  J>,Z)  so  that  the  sample  is 
parallel  to  the  electrode.  If  the  sample  were  not  parallel  to  the 
electrode  the  sensitivity  of  the  detector  circuit  would  be  markedly 
reduced.  The  can  is  then  put  into  place  and  the  chamber  evacuated. 

o 

The  sample  is  now  ready  for  the  measurements.  The  resonant  frequency 
is  found  by  connecting  the  input  of  the  driving  voltage  amplifier  to  a 
general  purpose  wide  band  oscillator  and  sweeping  the  frequency  through 
the  region  where  it  is  expected.  A  table  giving  the  characteristics  of 
the  first  flexural  vibration  modes  in  a  thin  circular  plate  is  given  in 
Appendix  I.  The  driving  oscillator  is  now  tuned  to  this  frequency  and 
switched  into  the  amplifier.  The  filter  is  then  switched  into  the  de¬ 
tection  circuit  and  the  tuned  amplifier  is  adjusted  for  maximum  gain. 
The  sample  is  now  being  driven  at  one  of  its  resonant ' frequencies  and 
these  vibrations  are  being  detected  and  displayed  on  the  oscilloscope 
along  with  the  driving  voltage. 


A  quick  check  at  this  point  will  indicate  whether  the  width 
of  the  resonance  peak  is  sufficient  to  allow  measurement,  or  not.  Since 
the  frequency  counter  will  only  allow  measurements  of  this  width  to 
within  0.2  H  it  is  necessary  that  it  be  several  H,_  in  order  to  get  an 
accurate  value  for  If  it  is  possible  to  make  this  measurement 
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then  by  using  the  vacuum  tube  voltmeter  and  frequency  counter  the  two 
frequencies  and  f£  at  which  the  amplitude  of  vibration  is  one  half 
that  at  resonance  are  measured  and  is  calculated  from, 


which  was  given  in  Chapter  II.  A  typical  resonance  curve  is  shown  in 
figure  3.6(a). 


When  the  resonance  peak  is  too  narrow  to  be  measured  with 


any  appreciable  accuracy,  which  occurs  'when  the  damping  is  low,  another 
method  must  be  used,  A  satisfactory  method  is  to  shut  off  the  driving 
voltage  and  measure  the  rate  of  decay  of  the  amplitude  of  vibrations 
at  resonance.  This  is  done  by  recording  the  detected  signal  on  the 
logarithmic  level  recorder  as  a  function  of  time.  The  trace  of  a  typi¬ 
cal  decay  is  shown  in  figure  3.6(b).  Since  the  decay  of  vibrations  is 
exponential  it  produces  a  straight  line  on  the  recorder.  This  line  is 
then  extrapolated  to  50  db  and  the  time  of  this  50  db  decay,  t^,  is 
then  determined  from  the  recorder  speed.  The  value  of  Q”  can  be  cal¬ 


culated  from  this  information  using: 


Figure  3 .6 


(a)  Resonance  curve  for  a  specimen  showing  characteristic  frequencies 

(b)  Typical  logarithmic  recording  of- the  decay  of  vibrations 
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which  was  derived  from  equation  (2.26)  in  chapter  II.  Satisfactory 
agreement  has  been  found  between  the  two  methods  when  it  was  possible 
to  use  both  at  the  same  time. 

It  should  be  emphasized  that  these  two  methods  are  comple¬ 
mentary.  When  the  damping  of  the  vibrations  is  high,  the  resonance 
peak  is  sufficiently  wide  to  yield,  an  accurate  value  of  Q"^.  As  the 
damping  gets  smaller,  the  peak  becomes  too  narrow  for  measurement, 
however  the  decay  time  becomes  long  enough  for  accurate  measurement, 
that  is  greater  than  0.5  seconds. 

The  sample  and  chamber  were  cooled,  from  room  temperature  to 
liquid  helium  temperature  and  measurements  were  made  at  intervals  of 
5°  to  10°  K  in  order  to  get  the  temperature  dependence  of  Q”~.  The 
cooling  was  achieved  by  placing  the  chamber  above  a  boiling  liquid  gas. 
Liquid  nitrogen  was  used  to  cool  to  77°  K  and  the  remainder  of  the 
temperature  range  was  covered  using  liquid  helium.  Thermal  contact 
between  the  sample  and  the  chamber  was  maintained  by  helium  exchange 
gas.  The  temperature  was  kept  constant  during  the  measurements  by 
pumping  out  this  exchange  gas.  At  the  beginning  and  end  of  each  me as- 
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urement , 
within  0 
helium. 


the 


.2° 


temperature  was  noted  and  it  was  found  that  the  drift  wa 
K  except  when  approaching  the  boiling  point  of  liquid 


o 


E.  Cryostat  Test 


In  order  to  estimate  the  damping  of  vibrations  due  to  the 
background,  that  is  the  mount  and  the  atmosphere  around  the  sample,  a 
trial  run  was  carried,  out  on  a  well  annealed  sample  of  tantalum.  This 
metal  was  chosen  because  of  its  very  low  intrinsic  damping  as  was  shown 
by  Bordoni  (1961).  The  results  of  these  measurements  are  shown  in  fig¬ 
ure  3.7.  As  expected  the  damping  decreases  with  temperature  and,  at 
about  15°  K,  Q“^  reaches  a  value  of  2.66  X  10“^.  The  damping  measured 
below  this  temperature  showed  considerable  amplitude  dependence,  which 
was  attributed,  to  damping  due  to  external  causes.  It  was  estimated 
that  external  contribution  to  damping  mads  up  3 0$  of  the  recorded  value, 

*)  n 

that  is  it  was  equivalent  to  a  Q  of  approximately  1  X  10“ ! ,  This  is 
less  than  Vjo  of  damping  measured  in  iron  at  its  minimum  value  (refer 
chapter  IV)  so  it  shall  be  neglected,  in  discussing  the  results. 
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Figure  3.7 


The  dependence  of  the  dissipation  coefficient,  Q“~,  on  temperature  in 
tantalum  (annealed,  in  vacuum  for  8  hrs.  at  800°  C) 
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CHAPTER  IV 


SXPER H' j3 NT AL  RESULTS  Al'JD 


DISCUSSION 


The  motivation  for  this  work,  as  it  was  pointed  out  earlier, 

c 

was  to  study  the  anomalous  rise  in  vibration  damping  observed,  at  low 
temperature,  in  iron.  Normally,  in  a  metal,  the  damping  of  vibrations 
slowly  decreases  with  temperature  having  its  lowest  value  at  0°  K 
(Mason,  1961).  The  first  st  ep  in  this  study  of  iron  was  to  measure 
the  temperature  dependence  of  Q“~.  The  graph  of  figure  4.1  is  a  plot 
of  the  results  obtained,  for  four  different  frequencies.  These  frequen 
cies  represent  the  fundamental  mode  f^  and  the  first  harmonic  fg  of 
flexural  vibrations  without  nodal  diameters  in  two  samples.  These  two 
samples  were  cut  from  the  same  plate  and  after  machining  they  were 
annealed  together  in  vacuum  for  eight  hours  at  400°  C.  At  the  end  of 
this  eight  hours,  the  temperature  of  the  samples  was  slowly  reduced  to 
room  temperature  by  cycling  it  around  a  decreasing  mean  value.  This 
method  of  annealing  was  found  to  give  the  lowest  room  temperature  damp 
ing  and  the  largest  low  temperature  anomally.  The  physical  dimensions 
of  all  the  samples  studied  in  the  work  are  given  in  Appendix  II. 

In  order  to  explain  this  anomalous  behavior,  its  dependence 
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on  the  amplitude  of  the  vibrations  must  also  be  determined.  The  dis¬ 
sipation  coefficient,  Q~^"  was  found  to  be  independent  of  the  amplitude 
of  the  vibrations  throughout  the  temperature  and  frequency  range  for 
all  amplitudes  used  in  this  study*  The  calculation  of  this  vibrational 
amplitude  range  was  described  by  Kuovo  (1961).  It  was  shown  in  chapter 
III  (equation  3»l)  that  the  force  applied  to  the  sample  was  t 


F  _  t°t  if 

r  —  p  «  « 

2d 


Taking  the  real  part  of  equation  (2.27)  we  find  that  the  amplitude  of 
the  vibrations  X  ,  is  given  by; 


C  J  VI2, 

O  0_ 

2  d^p&)  ^ 


Q 


at  the  resonance  condition.  In  order  to  apply  this  equation  to  the 
sample,  A  is  taken  as  the  area  of  the  electrode  and  p  as  the  density  of 
the  sample.  The  strain  amplitude  €  at  the  center  of  the  plate  will  then 
be  given  by: 

^  h  X 
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where  h  is  the  thickness  of  the  plate  and  a-^  is  the  radius  of  the  smal 
lest  nodal  circle  (see  Appendix  I). 

The  values  of  the  parameters  involved  were  within  the  fol¬ 
lowing  ranges: 

-3  -2 

5  *  10  cm  <:  d.  <:  5  x  13  cm, 

30  volts  <c  V  ^  100  volts, 

10  mm  a-^  ^  15  mm, 

103  <C  Q  ^  10-5 , 

60  KH  -<:  O  <  400  ifri 

z  o  z 

and , 

A  =  100  mm2, 

h  =  3 

f  ~  8g/c ir?. 

Taking  these  into  account,  it  is  found  that  the  range  of  the  experi¬ 
mental  strain  amplitude  is  approximately  given  by: 

-12  -Q 

5  X  10  <  £  <  3  x-  10  "  . 

•  • 

This  gives  a  maximum  value  to  X  of  the  order  of  a  few  Angstrom  units. 


'17 


It  can  ba  seen,  then,  that  the  vibration  amplitudes  being  studied  are 
of  the  same  order  of  magnitude  as  the  interatomic  distances.  Any 
attempt  to  use  vibrations  with  larger  amplitudes  than  that  mentioned 
above  caused  the  dissipation  to  increase  and  become  amplitude  dependent, 
This,  as  yet  unexplained,  behavior  is  thought  to  be  unrelated  to  the 
effect  being  studied  in  iron  since  it  occured  at  all  temperatures,  in¬ 
cluding  room  temperature  where  the  anomaly  is  found  to  be  negligible. 
Any  interpretation  of  the  anomalous  behavior  must  therefore  show  the  ■ 
amplitude  independence  of  the  dissipation,  at  least  for  low  amplitude 
vibratioiis. 

It  was  shown  in  Chapter  1  that  Q-^  is  a  measure  of  the  energy 
dissipated  by  vibrations  to  their  supporting  medium  and  it  was  pointed 
out  that  this  dissipation  may  occur  through  more  than  one  process. 


For  this  reason  it  will  be  considered  that  the  measured  values  of  Q' 
can  be  written  in  the  following  form: 


-1 


Cf1  =  +  AQ”1  , 


where  Q"^  is  the  damping  due  to  dissipation  of  energy  through  all  pro¬ 
cesses  except  that  which  causes  the  anomalous  rise  as  low  temperature. 
This  rise  is  noted  AQ~\  The  value  of  was  chosen  to  be  slowly 
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decreasing  and  linear  in  temperature.  This  is  illustrated  in  figure 
3.1  and  was  justified  at  the  beginning  of  this  chapter  by  what  is  found 
in  a  normal  metal.  An  example  of  this  behavior  is  shown  in  figure  3.^ 
which  gives  the  variation  of  with  temperature  in  nickel.  The  value 
of  AQ"^  can  now  be  found  by  simple  subtraction.  The  dependence  of  AQ  ^ 
on  frequency,  as  derived  from  figure  4.1,  is  shown  in  figure  4.2  for 
various  temperatures. 


In  his  paper,  Heller  (1961)  attempted  to  explain  AQ~1  in 
terms  of  the  magnetic  properties  of  iron,  more  specifically  the  magneto- 
strictive  properties.  It  was  his  contention  that  this  anomalous  behavior 
was  related,  the  the  AE  effect  described  by  Eozorth  (1951>  PP.  684).  If 
this  interpretation  is  to  be  valid,  then  AQ“^  should  be  dependent  on 
magnetic  field.  More  specifically,  the  damping  due  to  magnetostrietively 


induced  eddy  currents  should  go  to  zero  at  saturation  fields.  The 
data  illustrated  in  figure  4.3  show  no  appreciable  dependence  of  the 
damning  on  magnetic  field.  It  was  also  thought  advisable  to  measure 
the  damping  in  other  ferromagnetic  metals  in  order  to  compare  their  be¬ 
havior  with  that  of  iron.  The  results  for  nickel  are  shown  in  figure 
4.4.  There  is  no  indication  from  these  data  that  nickel  behaves  anoma¬ 
lously  like  iron  does. 


4  9 


The  variation 
frequency  for 


Figure  4.2 


of  the  anomalous  dissipation  coefficient,  AQ"*'1  ,  with 
several  temperatures  in  iron. 
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The  effect  of 
line  represen 


Figure  4.3 


a  magnetic  field  on  the  dissipation  in  iron.  The  solid 
ts  data  taken  with  no  magnetic  field. 
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Figure  4.4 

Damping  of  vibrations  in  nickel 

O  -  as  received  condition 

O  -  after  annealing  6  hours  at  450°  C  in  vacuum 
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The  conclusion  that  must  be  drawn  from  this  study  is  that 
the  behavior  of  iron,  with  respect  to  vibration  damping  at  low  tempera¬ 
tures  cannot  be  related  to  its  magnetostrictive  properties  as  was  pro¬ 
posed  by  Heller, 

“In  a  further  attempt  to  substantiate  this  statement,  measure¬ 
ments  of  Q~“  were  done  on  cobalt.  The  graph  of  figure  4.5  illustrates 
the  results  of  these  measurements  on  the  sample  as  received  from  the 
supplier  and  figure  4.6  shows  the  variation  of  the  resonant  frequency 
with  temperature.  The  large  peak  seen  at  about  120°  K  and  its  history 
dependence  is  believed  to  be  due  to  the  unstable  hexagonal  lattice  of 
the  cobalt.  It  was  pointed  out  by  Bozorth  (1961,  pp.  262)  that  the  low 
temperature  hexagonal  phase  is  only  approximately  ideal  close  packed. 
This  behavior  is  explained  by  the  fact  that  the  hexagonal  phase  is 
stable  only  at  temperatures  so  low  that  the  diffusion  of  atoms  is  ex¬ 
tremely  slow  and  they  cannot  form  a  perfect  lattice  after  they  have 
once  been  disrupted  by  the  transformation.  This  transformation  from 
cubic  to  hexagonal  structure  occurs  at  approximately  425°  C. 

The  belief  that  the  behavior  shown  in  figure  4.5  and  4.6  was 
due  to  instability  in  the  lattice  was  substantiated  by  measurements 
taken  after  the  sample  had  been  annealed  above  the  transition  temperatur 
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Figure  4.5 


The  damping  of  vibrations  in  cobalts 


(as  received  condition)  showing 


a  hysteresis  effect 


COBALT 


Figure  4.6 


Variation  of  the  resonant  frequency  of  the  cobalt  sample  with  tempera¬ 
ture  (as  received  condition) 
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at  550°  C  for  three  hours.  The  temperature  was  then  slowly  reduced 
through  the  transition  temperature  and  maintained  at  310°  C  for  another 
three  hours.  From  this  point  the  cooling  to  room  temperature  took  four 
hours.  The  results  of  these  measurements  are  shown  in  figure  4.7. 

They  do  not  show  the  hysteresis  effect  which  had  been  seen  previously. 
It  should,  also  be  noted  that  after  annealing  treatments  which  involved 
quick  cooling  through  the  transition  temperature  measurements  of 
exhibited  hysteresis  similar  to  figure  4.5* 

When  the  data  in  figure  4 .7  Are  compared  with  the  results  for 
iron,  it  is  seen  that  cobalt  does  not  behave,  at  low  temperature,  in 
any  way  analogous  to  iron.  This  final  result  ties  in  with  the  belief 
that  the  ferromagnetic  properties  of  iron  are  not  the  cause  of  its 
anomalous  behavior. 

The  significant  characteristic  of  the  temperature  dependence 
of  the  damping  vibrations  in  iron  is  its  resemblance  to  the  electrical 
conductivity.  A  look  at  figure  4.8  shows  that  AQ”^  is  directly  propor¬ 
tional  to  the  electrical  conductivity.  The  conductivity  curve  shown  is 
after  White  and  Woods  (1959)  using  0.265  jtfW  cm  as  the  value  of  the 
residual  resistivity.  Measurements  made  on  a  bar  of  iron,  cut  from  the 
same  sheet  as  the  samples  and  annealed  at  the  same  time,  showed  that 
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Figure  4,7 

Damping  (0)  and  resonant  frequency  (O)  of  cobalt  after  annealing  and 
slow  cooling  through  the  cubic  to  hexagonal  transition 
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Figure  4.8 


Comparison  of  the  dissipation  coefficient,  AQ~\  and  the  electrical 
conductivity  of  iron.  The  solid  line  represents  the  conductivity 
(after  White  and  Woods,  1959)  and.  the  open  circles  are  the  measured 
of  iQ"1. 
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this  is  not  an  unreasonable  value.  The  above  comparison  was  clone  for 
low  frequency  data  (28  KHz).  When  this  procedure  was  attempted  for 
the  data  taken  at  a  higher  frequency,  that  is  140  KHz,  there  was  a 
considerable  deviation  from  direct  proportionality  with  the  electrical 
conductivity.  This  is  illustrated  by  the  following: 


Q' 


-1 


a  /T=60°K 


.  i  ,J  ^ 


a  /T=5°K  , 


where  a  is  the  electrical  conductivity. 


The  fact  that  the  damping  of  vibrations  is  proportional  to 

the  electrical  conductivity,  at  least  at  low  frequencies,  indicates  that 

the  energy  of  the  vibrations  is  being  dissipated  to  the  conduction 

electrons.  It  is  now  necessary  to  find  the  mechanism  by  which  this 

dissipation  occurs.  The  involvement  of  the  conduction  electrons  hints 

at  the  electron-phonon  interaction  as  a  possible  mechanism,  however,  as 

it  was  pointed  out  in  a  paper  by  Verdini  (1962 ),  the  observed  dissipa- 

3 

tion  is  a  factor  on  the  order  of  10  greater  than  the  one  predicted  by 
the  free  electron  theory.  It  has  been  stated,  by  Mason  (1966,  Vol  IV  B, 
pp.303)  and  in  a  paper  by  Jones  and  Raynes  (1964)  that  the  free  electron 
theory  underestimates  the  dissipation  of  elastic  energy  due  to  the 
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electron-phonon  Interaction  in  metals  with  a  non -spherical  Fermi  sur¬ 
face  such  as  iron,  however,  two  orders  of  magnitude  was  the  maximum 
discrepancy  mentioned  and  this  was  for  the  high  frequency  range  (approx¬ 
imately  1  GHz  in  molybdenum).  The  electron- phonon  interaction  theory 
also  predicts  that  the  dissipation  should  be  directly  proportional  to 
the  conductivity  for  the  frequencies  being  studied  here  (Ehatia,  1967), 
This  then  does  not  explain  the  discrepancy  seen  at  140  KHz,  It  must  be 
concluded  from  the  above  discussion  that  another  mechanism  is  required 
to  explain  the  anelastic  behavior  of  iron  at  low  temperatures. 


In  order  to  determine  the  nature  of  this  mechanism  it  was 
necessary  to  know  the  influence  of  other  parameters  on  the  dissipation 
factor,  It  is  known  from  the  literature  that  the  effect  being 

studied  is  smaller  in  commercial  grade  iron  than  in  the  spec-pure  grade 
used  in  this  work  (Verdini,  1962),  and  that  the  anomalous  rise  is  les¬ 
sened  by  the  addition  of  hydrogen  to  the  iron  (Heller,  1961).  It  was 
also  shown  by  Brunner  (i960)  that  the  amount  of  damping  at  low  tempera¬ 
ture  is  dependent  on  the  dislocation  density  in  the  sample.  The  fact 
that  the  effect  is  smaller  in  hydrogen  charged  and  commercial  grade  iron 
can  be  attributed  to  the  smaller  electrical  conductivity  due  to  impurity 
scattering .  The  dependence  on  the  dislocation  density  cannot,  however, 


be  explained  in  this  wa, y 


An  increase  in  the  dislocation  density  is 
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expected  to  reduce  the  conductivity  and  thereby  reduce  the  low  tempera¬ 
ture  damping.  The  measurements  illustrated  in  figure  4.9  indicate  that 
this  is  not  the  case.  After  a  5 cp  compressional  deformation  of  a  sample, 
a  change  in  the  damping  was  observed.  This  change,  however,  instead  of 
being  uniquely  negative  as  expected,  was  positive  at  130  KHz  and  nega¬ 
tive  at  34  KHz.  This  result  is  thought  to  indicate  that  the  dislocation 
play  as  direct  a  role  in  the  damping  mechanism  as  do  the  conduction  elec 
trons. 

Theory  of  . Acoustic  Attenuation  in  Metals 

Due  to  Dislocations  Damped  by  Electrons 


A  process  by  which  dislocations,  induced  to  move  by  the  oscil¬ 
lating  stress  field  in  the  sample,  are  damped  by  the  conduction  elec¬ 
trons  is  thought  to  provide  a  reasonable  explanation  for  the  observed 
phenomena.  The  theory  governing  this  process  was  reviewed  by  Mason,  in 
his  book  (1966,  Voi  XV  B,  pp.  310)  and  will  be  presented  hers. 


It  has  been  proposed  by  Koehler  (1952)  that  many  aspects  of 
mechanical  damping  due  to  the  presence  of  dislocations  can  be  described 
by  a  model  which  assumes  that  the  dislocation  line  segments  act  like  a 
damped  vibrating  string,  pinned  at  both  ends.  In  this  model,  the  equa¬ 


tion  of  motion  of  the  displacement,  y,  of  a  dislocation  line  in  the 
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Figure  4.9 

The  effect  of  cold  working  on  the  dissipation  in  iron.  The  solid  lines 
represent  the  measured  values  of  Q""^  in  an  annealed  sample  and  the 
circles  represent  measurements  made  on  the  same  sample  after  a  'yjo  defor¬ 
mation. 

A  -  f  =  34  KHz 
B  -  f  =  130  KHz 
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presence  of  an  oscillating  stress  field  is  then: 


>2 
o  y 

C  ^  =  force  per  unit  length, 

bx 


where  A  is  the  effective  line  density  and  C  is  the  tension  in  the  dis¬ 
location  line.  The  force  given  on  the  right  hand  side  of  the  above 
equation  includes  the  driving  force  of  the  stress  field  and  the  damp¬ 
ing  force  which  is,  in  this  case,  due  to  the  conduction  electrons. 

The  force  on  the  dislocation  per  unit  length  is  taken  tc  be  ab  (Cot¬ 
trell,  1964,  pp.  6),  where  a  is  the  shear  stress  component  in  the  direc¬ 
tion  of  the  displacement  and  b  is  the  burgers  vector  of  the  dislocation. 
This  model  then  assumes  that  the  damping  force  is  of  the  viscous  type, 
that  is,  proportional  to  the  velocity  of  the  displacement  and  the  equa¬ 
tion  of  motion  is  then  written: 


+ 


a 

0  ay 

C  tvu=  ct  b 
ox  o 


cos  art, 


where  B  is  the  damping  force  per  unit  length,  per  unit  velocity. 


' 
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It  was  pointed  out  by  Oen  et  al.  (i960)  that  experiment ally 
all  types  of  dislocation  damping  show  a  broad  maximum  in  their  fre¬ 
quency  dependence.  It  will  be  assumed,  here  that  the  region  of  interest 
lies  near  to  this  maximum  and  that  under  these  conditions,  the  motion 
of  the  dislocations  is  damping  limited  thereby  making  the  inertia  term 
negligible.  Under  this  assumption  the  equation  of  motion  must  be  writ¬ 
ten: 

2 

B  ~  -  C  —  =  a  b  cos  at,  (4.1) 

H  ix2  ° 

with  the  boundary  conditions  for  a  dislocation  segment  of  length  1 

being: 

y(-i/2,t)  =  yW/2,t)  =  0.  (4.2) 

It  is  now  necessary  to  evaluate  the  constant  B  for  the  case  of  electron 
damping  of  dislocations. 

A  dislocation  is  surrounded  by  a  strain  field  and  as  it  moves 
through  the  metal  the  strain  at  any  point  changes  as  a  function  of 
time.  A  changing  strain  indicates  the  atoms  making  up  the  metal  are 
moving.  The  free  electrons  in  the  metal  will  then  be  imparted  with  this 
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sarae  motion  in  order  to  establish,  current  neutrality  since  the  particles 
of  the  lattice  are  positively  charged.  These  moving  electrons  can  now 
be  scattered  by  other  electrons  and  thus  the  energy  of  the  elastic  vi¬ 
brations  can  be  dissipated.  In  order  to  develop  a  theory  for  this  form  . 
of  damping,  the  free  electron  gas  is  considered  to  be  a  viscous  fluid  of 

c 

viscosity  tj. 

In  order  to  calculate  the  drag  coefficient,  B,  for  the  moving 
dislocation  segment  it  is  necessary  to  know  the  form  of  the  strain  field 
around  the  dislocation.  The  simplest  case  to  consider  is  the  screw  dis¬ 
location  because  of  the  symmetry  of  its  strain  field.  The  following 
treatment  is  after  Cottrell  (1953*  pp.  •  If  the  z  axis  of  a  cylin¬ 
drical  coordinate  system  is  chosen  to  be  in  the  direction  of  the  dislo¬ 
cation  line,  a  screw  dislocation  is  defined  such  that  the  displacement, 
(j,  of  the  lattice  from  its  ideal  position  has  the  form: 


Gl>  =  b(  0/2tt)  , 


(4.3) 


where  b  is  the  magnitude  of  the  Burgers  vector  and  9  is  zero  at  the  slip 
plane.  It  follows  from  this  that  the  only  stress  component  present  will 
be  the  shear,  Jqz>  given  by: 
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a  -  H  . 

Oz  r  bQ- 

Considering  equation  (4.3)  this  can  be  written: 

JQz  =  Pb/2lTr* 

and  the  strain  component  will  be  given  by: 

S9z  =  b/2rrr  , 


( 4. 4) 


(4.5) 


(4.6) 


All  other  strain  components  of  the  dislocation  field  are  zero. 

Wh en  the  above  strain  field  moves  through  the  viscous  electron 
gas,  energy  is  dissipated.  This  energy  loss  can  be  calculated  for  an 
elementary  cylindrical  volume  from: 

(T0zSgz)r  dr  d9  dz  =  (f,S0zSQz  +  flSQz)r  dr  dz  »  (4.8) 


where  Ta  is  the  stress  component  containing  an  elastic  term,  ,  and 
3z  “  uz 

an  anelastic  term,  rjSftr ,  due  to  the  viscous  drag  of  the  electrons. 

y  z 

It  is  now  necessary  to  calculate  S0z  for  all  points  around  a 
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dislocation  moving  with  a  velocity  u,  Considering  a  point  initially  at 
a  distance  r,  from  the  dislocation  line  and  at  an  angle  0  from  the  di¬ 
rection  of  motion,  the  rate  of  change  in  Sa  for  a  time  interval  dt  is 
given  by: 


dt  dt 


(4.9) 


where  r^  is  given  by; 


r2  ~  rl 


1  - 


(cos  0)u  dtl 


The  small  limit  of  the  above  expression  yield-:’ 


bu  cos  9 
2  nr^ 


(4.10) 


The  rate  of  energy  dissipation  per  unit  length  of  dislocation  can  now  be 
calculated  by  substituting  the  expression  for  SQ  given  by  equation 

V  Lx 

(4.10)  and  that  of  Sq^  given  by  equation. (4.6)  into  equation  (4.8)  and 
integrating  over  a  cylinder  of  unit  length.  This  gives: 
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2rr 


r 


v 


0 
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since  the  first  term  in  (4.8)  is  the  elastic  term  and  does  not  contribute 
to  the  integral.  In  order  to  evaluate  the  drag  coefficient  B,  the  energy 
dissipated  by  the  strain  field  must  be  equated  to  the  work  done  by  the 
force  per  unit  length,  F,  applied  to  the  dislocation.  In  the  coordinate 
systexU  used  above, 


F  =  T13b  , 


(4.12) 


and  the  velocity  u,  attained  by  the  dislocation,  is  given  by: 


B 


(4.13) 


It  is  sesn  froin  equations  (4.11),  (4.12)  and  (4.13)  that 


(4.14) 


hence : 
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(4.1 5) 

The  quantity  aQ  appearing  in  equation  (4.1 5)  and  in  equation  (4.11)  is 
the  smallest  radius  around  the  dislocation  line  for  which  the  strain 

o 

field  is  elastic.  Inside  this  radius  the  motion  of  the  dislocation  is 
discontinuous  and  the  damping  due  to  the  viscosity  of  the  conduction 
electrons  is  assumed  to  be  zero.  The  dissipation  of  energy  via  acousti 
radiation  due  to  the  discontinuous  motion  of  the  atoms  within  this 
radius  is  neglected  in  the  theory.  A  value  of  3/4  b  has  been  given 
for  a^  by  Cottrell  (1953>  pp.  39) •  A  completely  analogous  treatment  of 
edge  dislocations  (xiason  1966,  Vol  I V  A,  pp.  313)  leads  to  a  similar 
expression  for  the  damping  constant  of  edge  dislocations .  In  this  case 


0.75  b2')  . 

8n(l-ap)2a2 


(4.16) 


where  <i  is  the  Poisson’s  ratio. 

P 

The  viscosity  q,  of  the  electron  gas,  is  evaluated  by  assum¬ 
ing  it  to  be  analogous  to  a  classical  viscous  fluid.  It  is  well  known 
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for  this  case  that, 


N  ‘miV 


3 


rms 

—  j 


(4.i?) 


where  N*  is  the  free  electron  density,  m  is  the  mass  of  the  electrons, 

t  is  the  mean  free  path  between  collisions  and  V  is  the  root  mean 

c  rms 

square  velocity.  In  the  case  of  a  metal  with  a  spherical  Fermi  surface, 
Vrms  can  4>e  replaced  by  3/5  of  the  Fermi  velocity  v^: 


n  = 


N'mv  4 
F 

5 


(4.18) 


When  the  mean  free  path  and  the  Fermi  velocity  are  expressed  in  terms  of 
the  electrical  conductivity,  a,  and  the  free  electron  density  (Kittel, 
1956,  chapter  10),  it  is  found  that: 


2,  2 

1i  (3n  iY 


'3/2 

y  j 


(4.19) 


where  e  is  the  electronic  charge  and  ¥\  is  Planck's  constant  divided  by 

2tt. 
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The  constant  C  in  equation  (4.1),  that  is  the  effective  ten¬ 
sion  in  the  moving  dislocation  line,  remains  to  be  evaluated.  It  has 
been  stated  by  Cottrell  (1964,  pp.  53)  that  an  acceptable  form  for  the 
tension  is  given  by  the  equation: 


(4.20) 


where  p  is  the  shear  modulus  of  the  metal. 

The  equation  of  motion  of  the  pinned  dislocation  segments, 


B 


'at  “ 


C 


iz 

6x2 


jQb  cos  a't  , 


(4.1) 


is  now  completely  known  since  the  constants  B  and  C  can  be  evaluated 
using  equations  (4.15),  (4.16),  (4.19)  and  (4.20).  In  their  paper,  Oen 
et  al.  (I960)  gave  the  steady  state  solution  of  equation  (4.1)  as: 


Q  b 
o 


r 


y(x,t)  =  __  ^ 


sin  cot 


V 


1  - 


cosh  kx  cos  kx  cosh(™)cos(~~i) 
_ _ _ _ 2 _ 2 

cosh^(Jsi)  -  sin^(iil) 

2  2 
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sinh  kx  sin  kx  sinh(™)sin(~^) 

,2 /kH \  .  2/kP  \ 

cosh  )  -  sin  (~~- ) 


+  cos  cot 


cosh  kx  cos  kx  sinh(^r)sin(^r) 


1 2  fk l \  .  2 ( k  £ \ 

cosn  {-—)  -  sm  (™) 


sinh  kx  sin  kx  cosh (■—-) cos {-—) 

£  C+ 

,  2/kP. \  •  2 / k^.\ 

cosh  (—)  -  sm  (~~) 


(4.21) 


2 

where  k  =  gjB/2C.  Expressed  as  the  sum  of  harmonics,  the  displacement  is, 


4a  b  °°  (  sin 

y(x,t)  =  y~j 


rr 


n  =  0  2n  + 1 


_  2 

a  (x+£/2)  /  a  C  cos  ot  u)B  sin  ^t 

L  n  J  /  n 

\ 


+ 


+u)2B^  a^C2  +co2b^ 


n 


(4.22) 


where  a^  =  (2n  +  l)rr/i,  Using  equation  (2.21)  the  dissipation  factor, 
Q"^,  is  calculated  from: 


OO 


,-l 


E(2)W(U)  dX 


1  o 


n 


azk 


(4.23) 
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where  W(.d)  dJ2-  is  the  number  of  dislocation  segments  per  unit  volume, 
with  the  distance  between  their  pinning  points  being  between  Z  and 
£+  dfl.  The  energy  loss  per  cycle  for  one  dislocation  segment  of  length 
i  is  given  by  E (A)  which  is  calculated  from: 


i/2 


2tt 


E(i)  = 


dx 


-i/2 


*  dy(x,t) 

a  b  cos  fut - dt 

°  dt 


(4.24) 


0 


Assuming  that,  the  distance  between  pinning  points  on  the  dislocations 
is  distributed  exponentially  about  a  mean  value,  .2,  (Koehler,  1952), 
that  is 


L  -i/l° 
NU)  =  —z  e  , 


(4.25) 


where  L  is  the  total  length  of  dislocations  per  unit  volume,  then  the 
value  of  Q"1  is  given  by 


/ sinh  kZQx  +  sin  kiQx\ 
(cosh  k l  x  +  cos  k 2-Qxj 


e”x  dx  c  (4.26) 
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In  order  to  account  for  the  fact  that  the  applied  stress  is 
not  always  of  the  form  and  direction  given  in  equation  (4.12),  it  is 
necessary  to  introduce  an  orientation  factor  R  (-/Hers  and  Thompson,  1961) 
so  that  equation  (4.26)  now  reads: 


Q' 


-1 


OO 


im: 


(keo)- 


0  L 


/sinh  ki  x  +  sin  k2  x  \ 

“V  -  / - rn - _2L  j 

1  cosh  kSL  x  -1-  cos  ki  x  / 


e-x  dx 


(4.27) 


The  range  of  R  is  given  by  the  above  authors  as  0.08  to  0.2.  The  dis¬ 
sipation  coefficient,  Q”"^,  as  given  by  equation  (4.27)  has  the  amplitude 
independence  that  is  required  by  the  experimental  results.  The  asymptotic 
behavior  of  the  above  equation,  for  low  frequencies,  is  given  by: 


Q'1  =  4  MR lZ  ~  ,  «  uj 


(4.28) 


and  for  high  frequencies  by: 


-1  2  ^4 

«  =  «“0  -£  - 


(a)  »  cJ  ) 
v  o' 


(4.29) 


where  ojQ  is  given  by: 


- 
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(4.30) 


From,  equations  (4.15),  (4.16)  and  (4.19)  it  is  seen  that  co 
must  be  inversely  porportional  to  the  electrical  conductivity.  The 
dissipation  of  elastic  energy  at  low  frequencies,  as  indicated  in  equa¬ 
tion  (4.28),  is  directly  proportional  to  the  electrical  conductivity. 

This  is  in  agreement  with  the  experiinental  results  which  are  illustrated 
in  figure  4.8.  The  deviation  from  this  direct  proportionality,  seen  at 
higher  frequencies,  is  also  expected  from  the  theory.  The  integral  of 
equation  (3.26)  has  been  calculated  numerically  (Oen  et  al. ,  i960)  for 
the  frequency  range  around  t'o.  The  result  of  this  calculation  is  the 
universal  curve  shown  in  figure  4.10.  It  was  possible  to  fit  the  experi¬ 
mental  results  illustrated  in  figure  4.2  to  this  curve.  The  best  agree¬ 
ment  between  theory  and  experiment,  as  illustrated  in  figure  4.11,  was 

o  2 

achieved  by  using  a  value  of  I.83  x  lO--5  for  the  product  and  a  value 

of  O  which  was  temperature  dependent.  Figure  4.12  is  a  plot  of  a>“ 
versus  temperature  as  used  in  this  fit.  It  is  now  necessary  to  compare 
these  parameters  with  their  accepted  values  in  iron.  According  to  Van 
Euren  (196l,  pp.  328),  a  reasonable  value  for  the  pumber  of  dislocations 

7  -2 

cutting  a  surface  of  unit  area  in  a  well  annealled  metal  is  10  cm  . 
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Figure  4.10 


The  frequency  of  Q“^  for  dislocation  damping  by  electrons  (ai 
et  al.,  I960) 
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Figure  4.11 


The  normalized  frequency  dependence  of  Q  in  iron 


Observed  results 


Theoretical  curve  after  Oen  et  al.  (i960) 


I 


i 

|i 


I 


I, 

i 


77 


Figure  4.12 


temperature  as  used  to  plot  observed  results  in  figure  4.11 


-electrical  conductivity  after  White  and  Woods  (1959) 
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This  is  equivalent  to  N,  the  total  length  of  dislocations  per  unit 
volume.  If  a  value  of  0.1  is  assumed  for  the  orientation  factor  R,  since 


=  1.83  X  10“3  , 

it  is  found  that  SL  must  have  a  value  of  4.3  X  lO”^  cm.  This  is  con¬ 
sidered  to  be  a  reasonable  average  length  for  dislocations  (Mason,  1966, 
Vol  IV  A,  pp.  320). 

In  order  to  compare  the  experimental  values  of  0)  with  theory, 
it  is  necessary  to  develop  further  equation  (4.30)  which  reads: 


U) 

o 


_2C_ 

Bl2  ' 


From  equation  (4.20)  it  is  seen  thatj 


(4.31) 


If  it  is  assumed  that  3  is  given  by  its  expression  in  equation  (4.16) 
for  edge  dislocations,  equation  (4.31)  can  be  written: 
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8tt^(1  -  a/a20 
0.75  J?|t| 


Using  equation  (4.19)  it  can  be  seen  that  ; 


(4.32) 


60  0  =  17.7 


/  \2  2 

^  v 

M»2/3g2f12J 
16  o 


(4.33) 


Assuming  that  the  effective  density  of  free  electron,  N‘ ,  is  independent 
of  temperature  and  since  the  elastic  modulii  are  very  nearly  independent 
of  temperature  (Bhatia,  1967),  it  is  found  that  (0  should  be  inversely 
proportional  to  the  electrical  conductivity  a.  The  experimental  values 
of  6 were  compared  with  the  electrical  conductivity  in  figure  4.2  and 
are  seen  to  be  in  very  good  agreement  with  this  finding.  Although  the 
temperature  dependence  of  u>  found  experimentally  is  in  agreement  with 
equation  (4.33) »  there  is  a  discrepancy  in  the  magnitude  of  this  param¬ 
eter.  In  order  to  evaluate  CO  from  equation  (4.33) »  a  value  of  0.2  free 
electrons  per  atom  was  used  to  arrive  at  N*  (Kott  and  Jones,  1958,  pp.  222). 
The  dislocations  were  assumed  to  all  be  of  the  unit  type  giving  aQ  a 
value  of  approximately  10”^  cm  (aQ  «  3/4  b).  This  led  to: 
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theory  =  9  *  (jaP-  cm  sec)  }  (4.34) 

From  figure  (4.12)  it  is  seen  that: 

(tJoCJ)exp.  =  2-45  x  107  (pift  cm  sec)"1  (4.35) 

It  is  believed  that  the  discrepancy  between  equations  (4.34) 
and  (4.35)  can  be  explained  through  either  of  two  weaknesses  in  the 
theory.  These  weaknesses  lie  in  the  calculation  of  the  damping  constant 
B  and  in  the  assumed,  distribution  of  dislocation  lengths  of  equation 
(4.25).  liquations  (4.15)  and  (4.16)  are  two  expressions  for  the  con¬ 
stant  B,  one  for  screw  dislocations  and  the  other  for  edge  dislocations. 
Although  both  of  these  expressions  give  the  required  temperature  depend¬ 
ence  of  o)  when  substituted  into  equation  (4.31),  that  of  equation  (4.16), 
giving  the  damping  constant  of  edge  dislocations  was  used  to  arrive  at 
the  value  of  in  equation  (4.34).,  the  reason  being  that  better  agree¬ 
ment  was  obtained  in  this  way  between  expr  essions  (4.34)  and  (4.35)*  It 
is  possible  that  a  third  type  of  dislocation  is  responsible  for  the  ob¬ 
served  behavior.  In  applying  the  damping  constant  of  edge  dislocations 
to  equation  (4.31)  it  was  assumed  that  all  the  dislocations  were  of  the 
unit  type.  This  was  not  done  in  consideration  of  the  properties  of  iron 


■ 
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but  only  to  facilitate  calculations.  This  assumption  fixed  the  value 

-8 

of  a^  at  10”  cm.  It  is  very  probable  that  the  Burgers  vectors  of  the 
dislocations  in  iron  have  a  range  of  values  and  it  is  possible  that  be¬ 
cause  of  this  the  parameter  aQ  has  a  smaller  effective  value  than  the 

one  assumed  above.  This  would  lead  to  a  smaller  theoretical  value  for  u)  . 

o 

With  reference  to  the  distribution  of  the  distance  between 
pinning  points  assumed,  in  equation  (4.25),  it  has  been  shown  by  Gen  et  al. 
(I960)  that  if  one  assumes  a  delta  function  distribution  for  this  length 
that  the  peak  shoxm  in  figure  4.10  is  shifted  to  a  larger  value  of  ayu^. 
The  general  characteristics  of  the  temperature  and  frequency  dependence 
of  Q“  ,  however,  remain  much  the  same.  In  this  case,  a  fit  of  the  ex¬ 
perimental  data  like  that  of  figure  4.11  would  give  an  even  smaller  ex¬ 
perimental  value  to  c J  .  The  delta  function  distribution  of  the  distance 
between  pinning  point  must  be  discounted  because  of  their  random  nature, 
however,  a  distribution  which  would  produce  a  shift  of  the  peak  in  the 
opposite  direction  cannot  be  ignored.  It  is  well  known  that  the  dis¬ 
location  structure  of  iron  is  peculiar  to  this  metal  as  evidenced  by  the 
fact  that  no  Bordoni  peak  has  yet  been  observed  in  it.  This  may  also 
explain  the  fact  that,  in  metals  which  have  similar  crystal  and  electronic 
structures  to  that  of  iron  such  as  niolium  and  tantalum,  there  have  been 
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no  observations  of  this  type  of  damping.  In  summary,  it  is  thought  the 
fact  that  6-)q,  as  measured  in  iron,  is  much  smaller  than  expected  can  be 
explained,  by  a  peculiar  dislocation  structure  in  iron. 

This  small  value  of  provides  an  explanation  for  measure¬ 
ments  done  in  the  megacycle  frequency  range  by  Verdini  (private  communi¬ 
cation).  During  the  course  of  this  work  the  type  of  damping  which  occurs 
at  lower  frequency  was  not  observed.  For  this  frequency  range  it  is 
expected  from  equation  (4.29)  that  Q~“  will  be  directly  proportional  to 
the  electrical  resistivity.  The  behavior  would  be  almost  indistinguish¬ 
able  from  the  background  damping  due  to  dislocation-phonon  and  dislocation' 
dislocation  interactions. 
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CHAPTER  V 


CONCLUSIONS  AND  SUGGESTION'S  FOR  FURTHER  WORK 


The  anomalous  behavior  of  vibration  damping  in  iron,  at  low 
temperatures,  has  been  studied  experimentally.  The  proportionality  of 
the  dissipation  coefficient,  Q~^,  to  the  electrical  conductivity  and 
its  dependence  on  the  density  of  dislocations  led  to  an  interpretation 
of  this  effect  as  being  due  to  the  damping  of  the  moving  dislocation 
lines  in  the  sample  by  the  conduction  electrons.  The  theory  of  this 
damping  mechanism  is  presented  in  the  text  and  good  agreement  is  shown 
to  exist  between  theory  and  experiment  in  the  form  of  the  temperature 
and  frequency  dependence  of  Q”^.  The  theory  presented  indicates  that  a 
maximum  value  of  should  occur  at  a  certain  frequency.  The  value  of 
this  frequency  determined  from  the  experimental  results  is  much  smaller 
than  that  predicted  by  the  theory.  It  is  thought  that  this  discrepancy 
can  be  explained  by  the  fact  that  the  theory  does  not  take  into  account 
the  possible  existence  of  a  dislocation  structure  which  is  peculiar  to 
iron. 

In  order  to  clarify  the  existing  situation,  a  series  of  ex¬ 
periments  have  been  planned  for  the  future. 
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(1)  A  study  of  the  effect  of  y  irradiation  on  the  dissips.tion  coefficient 
Q  ^  would  provide  a  controlled  method  of  studying  the  effect  of  changes 

in  the  dislocation  structure  of  iron.  In  complement  to  this  a  systematic 
study  of  the  effect  of  various  annealing  treatments  and  cold-working 
could  provide  new  experimental  information,  modification  of  the  existing 
equipment  is  being  considered  in  order  to  allow  measurements  of  both 
and  the  electrical  resistivity  in  the  same  sample.  Any  changes  in  the 
resistivity  due  to  the  above  treatments  could  then  be  taken  into  account 
when  studying  the  dislocation  structure.  This  is  necessary  since  both  of 
these  properties  effect  the  value  of 

(2)  Since  a  sample  vibrating  in  its  flexural  mode  experiences  both  shear 
and  longitudinal  stress,  it  is  not  known  whether  vibrations  which  contain 
only  pure  shear  or  pure  longitudinal  stresses  are  damped  in  the  same  way. 
For  this  reason,  measurement  of  for  extensional  and  torsional  vibra¬ 
tions  in  a  bar  shaped  sample  of  iron  are  planned  for  a  wide  range  of 
frequencies. 

During  the  course  of  the  present  research  the  dissipation 
factor,  Q”1,  was  measured  in  cobalt  and  nickel.  These  two  metals  did 
not  exhibit  the  behavior  seen  in  iron.  The  measurements  on  cobalt,  how¬ 
ever,  did  show  evidence  of  the  unstable  low  temperature  hexagonal  phase 
of  this  metal. 
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APPENDIX  I 


FARAI4ET.ERS  FOP 


FLEXURAL  VIBRATIONS  IN  THIN  PLATES* 


For  a  circular  plate  of  radius  a  and  thickness  h,  the  fre¬ 
quency  of  the  normal  flexural  modes  of  vibration  are  given  by 


^m,n  ~  ccm,n 


E  |1/Z 
\p(l-o2)/ 


(1-1) 


where  E  is  Young’s  modulus,  a  is  the  Poisson’s  ratio  and  o  is  the  density 

P 

of  the  material.  The  constants,  a  are  the  eigenvalues  of  the  modes 

m,n  b 

with  m  nodal  circles  and  n  nodal  diameters.  The  values  of  a  _  are 

m,  n 

given  in  the  following  table: 


III 

n  =  0 

n  =  1 

n  =  2 

n  =  3 

0 

- 

0.241 

0.562 

1 

0.417 

0.943 

1.619 

2.431 

2 

1.770 

2.750 

.  • 

_ 

*  After  Nuovo  (1961) 
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Since  the  force  of  the  electrode  on  the  sample  was  applied 
almost  uniformly  over  a  circular  area,  only  modes  with  this  symmetry 
were  used  for  measurements,  that  is  modes  with  no  nodal  diameters.  The 
radius  of  the  nodal  diameters  for  the  first  two  such  modes  is  given  by: 


r-^  q  =  0.673a; 


f 


2,0 


‘2,0 


0.392a 

0.842a 


(1-3) 


When  running  in  the  second  mode  the  mounting  pins  were  placed  on  the 
nodal  circle  of  radius  r"  since  it  is  larger  than  r '  and  affords 

t  )\J 

better  support  for  the  sample. 


APPENDIX  II 


PxiYSICAL  CHARACTERISTICS  OF  SAPLBS 

IRON  (Johnson  iiatthey  spectr ©graphically  standardized) 
A.  Spectral  .Analysis 


Element 

Quantity  present  (ppm) 

Silicon 

8 

Manganese 

3 

Nickel 

3 

Copper 

2 

Magnesium 

1 

Silver 

less  than  1 

B.  (1)  Figure  4.1 


Sample  #1:  diameter  - 

30  mm 

thickness 

-  30  mm 

Sample  #2 :  diameter  - 

25  mm 

thickness 

-  3  mm 

(2)  Figures  4.3  and  4.9 

Before  compression:  diameter  -  25  mm 


thickness  -  2.6  mm 
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After  compression:  diameter  -  25  mm 

(using  force  of  thickness  -  2.5  mm  (irregular) 

25,000  lbs.) 

NICKEL  (source  unknown ) 

Figure  4.4:  diameter  -  36  mm 
thickness  -  5  mm 

COBALT  (Johnson  Latthey  spectro graphically  standardized) 

A.  Spectral  Analysis 

Impurity  content  found  to  be  very  similar  to  that  of  iron,  that  is 
less  than  20  parts  per  million. 

B.  Figure  4.5,  4.6  and  4.7:  diameter  -  30  mm 

thickness  -  5  mm 


